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We study the effects of electron correlation on transport through an interacting region con- 
nected to multi-mode leads based on the perturbation expansion with respect to the inter- 
electron interaction. At zero temperature the conductance defined in the Kubo formalism 
can be written in terms of a single-particle Green's function at the Fermi energy, and it can 
be mapped onto a transmission coefficient of the free quasiparticles described by an effective 
Hamiltonian. We apply this formulation to a two-dimensional Hubbard model of finite size 
connected to two noninteracting leads. We calculate the conductance in the electron-hole sym- 
metric case using the order U 2 self-energy. The conductance shows several maximums in the 
U dependence in some parameter regions of t y /t x , where t x (t y ) is the hopping matrix element 
in the x- (y-) directions. This is caused by the resonance occurring in some of the subbands, 
and is related with the U dependence of the eigenvalues of the effective Hamiltonian. 
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§1. Introduction 

Low-dimensional electron systems have been one of the current interests in the fields of the 
condensed matter physics and materials science. For instance, in some of the organic conductors, 
the electron correlation has been considered to be important to understand the physical properties.^ 
The Kondo effect in quantum dotsSH& has also been studied intensively from theoreticalS!** and 
experimentalBili) sides. When the average number of the electrons in a quantum dot is odd, 
the perfect transmission due to the Kondo resonance situated at the Fermi energy occurs at low 
temperatures. Recently, artificial molecules which are realized by arranging two or more quantum 
dots have also been studied0ini) Theoretically, the crossover from the high-temperature Coulomb- 
blockade to low-temperature Fermi-liquid behaviors of the quantum dots has been studied using 
advanced numerical methods such as the numerical renormalization groupQi) and quantum Monte 
Carlo methodsSii© 

In a previous work, one of the authors has studied the conductance of small interacting systems 



connected to two single-mode leads Jll'lll' and has calculated the conductance of a Hubbard chain 
of finite size N using the order U 2 self-energy. The results obtained in the electron-hole symmetric 
case depend strongly on whether N is even or odd. In the even cases, the conductance decreases 
with increasing N showing a tendency toward a Mott-Hubbard insulator. On the other hand, in 
the odd cases the perfect transmission due to the Kondo resonance occurs© 

The purpose of the present work is to generalize the formulation to the multi-mode systems where 
the interacting system is connceted to noninteracting leads with a number of channels. As in the 
single-mode case, at T = the contributions of the vertex corrections on the dc conductance g 
vanish. Then the conductance is determined by the value of the single-particle Green's function at 
the Fermi energy u) = 0, and thus g can be written in terms of the transmission coefficient of the 
free quasiparticles described by an effective Hamiltonian also in the multi-mode case. We apply 
the formulation to a two-dimensional Hubbard model of finite size consisting of Nq = N x M sites, 
where N and M are the size in the x- and y-direction, respectively. Two noninteracting leads of 
M channels are attached to this cluster along the x-direction. This system may be considered as 
a model for two-dimensional materials such as an array of quantum dots and a carbon nanotube. 
For instance, for M = 2 it may also be considered as a model for a ladder system of nanometer 
size. Since there is no translational symmetry in the systems we are considering, the self-energy has 
Nq x Nq matrix elements. We calculate all these elements within the second order perturbation 
expansion in U in the electron-hole symmetric case, and determine the parameters of the effective 
Hamiltonian up to terms of order U 2 . The results of the conductance show maximums at finite U 
for some values of t y /t x , where t x (t y ) is the nearest-neighbor transfer integral in the x-direction 
(y-direction). This is caused by resonances occurring in some of the subbands. These behaviors can 
also be understood through the U dependence of the eigenvalues of the effective Hamiltonian. We 
note that the second-order perturbation theory has been applied by a number of groups to study 
transport through the single Anderson impurit and systems consisting of a number of 

sites.S'S-* Also, the band calculations has been applied to obtain the conductance of nanometer 
systems.il' Our approach can, in principle, be extended to realistic systems by evaluating the order 
U 2 self-energy using the orbitals which are obtained with the band calculations. 

In §2, we describe the outline of the general formulation and introduce the effective Hamiltonian. 
In §3, we apply the method to a two-dimensional Hubbard model, and present the results of the 
conductance. Summary is given in §4. In Appendix, we provide the derivation of the expressions 
of the conductance and the total charge displacement. 

§2. Formulation 

In this section, we give the outline of a general formulation of the conductance through inter- 
acting systems connected to Fermi-liquid reservoirs with a number of channels. Our formulation 
is applicable for various systems which have the time reversal symmetry, i.e., eqs. (|2.9|) and (2.11) 
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hold for a wide range of the systems described by the Hamiltonian eq. (2.1). We provide the details 
of the derivations in Appendix, and demonstrate the application to a Hubbard model in the next 
section. 

We start with a system which consists of three regions as illustrated in Fig. 1; a finite interacting 
region (C) situated in the middle, and two noninteracting reservoirs on the left (L) and right 
(R). The central region contains iVc resonant levels, and the inter-electron interaction Uj 4 j 3; j 2 j 1 is 
switched on in the central region. Each of the two reservoirs is infinitely large and has a continuous 
energy spectrum. The central region and two reservoirs are connected by Ml and Mr channels 
which are described by the mixing matrix elements v L m and v R m , respectively. The complete 
Hamiltonian is given by 

Ti = Hh + Ti-R + T~ic + ^C* + ^mix , (2.1) 
ijeL (J 
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Here cj CT (c jo .) creates (destroys) an electron with spin a at site j, and fi is the chemical potential. 
t\p tfj, and i^j- are the intra-region hopping matrix elements in each of the regions L, R, and C, 
respectively. The labels 1, 2, . . ., Afc are assigned to the sites in the central region. In eq. (2.6), T m 
(£ m ) is the label assigned to the m-th site at the sample side (lead side) of the interface on the left. 
A/" m (1Z m ) is the label assigned to the m-th site at the sample side (lead side) of the interface on 
the right. Note that the number of the channels is less than the number of the interacting sites, i.e, 
Ml < Nq an d Mr < Nq- For instance, in the case of Ml = Mr = Nq, all the sites in the central 
region are connected to both of the reservoirs, i.e., X m — ■N'm for m — 1; 2, . . . , Nq. We assume 
that all the hopping matrix elements are real, and the interaction has the time reversal symmetry: 
^43;2i is real and = Uu-u = £/i2;34 = £/42;3i = t r i3 ; 24- We will be using units U = 1 unless 

otherwise noted. 
The single-particle Green's function is defined by 

GjAien) = - j^&r (T T c ja (r) c], a (0)) e i£ ^. (2.7) 
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Here (3 = 1/T, e n = (2n + 1)tt/(3, Cj a (r) = e T ™Cj a e rW , and (• • •) denotes the thermal average 
TV e~P H --- /Tre-P n . The spin index has been omitted from the left-hand side of eq. ( |2.7| ) 
assuming the expectation value to be independent of whether spin is up or down. Since the 
interaction is switched on only for the electrons in the central region, the Dyson equation can be 
written as 

Gij(z) = G%{z) + G iA z ) G fj (z) . (2.8) 

i'j'ec 

Here G^{z) is the unperturbed Green's function corresponding to the noninteracting Hamiltonian 
7i° = TCl + "Hr + + Ttmix- The summations with respect to i' and j' run over the sites in 
the central region, and Y^iiji(z) is the self-energy due to the interaction TCq 1 . Because of the time 
reversal symmetry of Ti, these functions are symmetric against the interchange of the site indices: 
Gij(z) = Gji(z) and Sjj(z) = Sjj(z). Note that at T = the single-particle excitation at the 
Fermi energy z = i0 + does not decay, i.e., ImS+.(0) = oH> In what follows, we will treat 
complex variable, and use the symbol + (— ) in the superscript as a label for the retarded (advanced) 
function: Ej(e) = E# (e ± i0+). 

As in the single-mode case, the dc conductance at T = can be expressed in terms of the 
single-particle Green's function at u = [see Appendix [~A| ; 

9 = ^-Tr [4r R (0)G^(0)r L (0) G?^(0)] . (2.9) 

Here -Tl(w) and -Tr(w) are Ml x Ml and Mr x Mr matrices, respectively. These two matrices 
are caused by the coupling with the left (L) and right (R) leads, and the elements are given 
by T a . mm ,(uj) = -Im v ajm F+ mm , (uj) v aiTn > with F a _ mm /(z) being the Green's functions at the 
interface of the isolated lead (a = L,R). In eq. ( |2.9| ), Gj^ x and Gjj^ are Mr x Ml and Ml x Mr 
matrices the elements of which are given by Gj^ iIm and G^^y^, respectively. 

Another quantity which can be related to the scattering coefficients is the displacement of the 
total charge defined byEl'0 

Aiv tot = EE(4o 
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40 - (clc ia ) R . (2.10) 



Here (■ ■ -) L and (• • -) R denote the ground-state average of isolated leads described by Ti.^ and 
7Yr, respectively. At T = 0, following the derivation of the Friedel sum rule by Langer and 
Ambegaokar,© AA^ot can be written in terms of a (Ml + Mr) x (Ml + Mr) matrix S, 

AiVtot = - log[det5] , (2.11) 

7T1 
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Here Gjj and G~j^j^ are Ml x Ml and Mr x Mr matrices, and the elements are given by Gj iXm 
and Gj^j^ m , respectively. The outline of the Friedel sum rule in the single-mode case is given in 
Appendix [^. The extension to multi-mode case eq. (2.11) is straightforward because in the case 
of the Friedel sum rule we do not have to take into account the contributions vertex corrections. 
For convenience, in eq. (2.12) we have assumed that there is no common element in the sets X 
and J\f. This is not an important assumption, and the expression without the assumption can be 
obtained in the similar way. Note that eq. (2.11) is also written in terms of the Green's function 
at T = and u = 0. Thus, due to the property ImS^(O) = 0, the values of g and AiV tot at T = 
can be expressed in terms of the transmission and reflection coefficients defined with respect to a 
one-particle Hamiltonian@ ) 

Wqp = T~(-L + T~(-R + %C + T^mix • (2-13) 

Here 1~Cq = YlijeC,a ( ~Wj ~ M<% ) c la c ja w ^ n ~~ Wj = + ^ e ^ij(0)- This effective Hamiltonian 
describes free quasiparticles of the local Fermi-liquid. As shown in the next section, the eigenvalues 
of the partial Hamiltonian T~Cq have important information about the ground-state properties. 

§3. Application to a Two-dimensional Hubbard Model 

In this section, we apply the formulation to a two-dimensional Hubbard model connected to reser- 
voirs and calculate the conductance using the order U 2 self-energy in the electron-hole symmetric 
case. 

3.1 Model 

The schematic picture of the model is illustrated in Fig. 2. The system size in the y-direction, M, 
is finite. In the x-direction, the central region consists of ./V columns, and two noninteracting leads 
are connected at x = 1 and x = N. Thus, the total number of the interacting sites is Nq = N x M. 
The parameters of the Hamiltonian eq. (2.1) are specified as follows. The mixing matrix element 
is assumed to be uniform in the y-direction: u Lm = d l and v Km = v R for m = 1, 2, . . . , M. We 
assume that the off diagonal part of tfj describes the nearest-neighbor hopping; t x and t y in the 
x- and y-directions, respectively. Along the y-direction, we use the periodic boundary condition. 
Furthermore, we assume Uj 4 j 3] j 2 j 1 to be an onsite repulsion U, and concentrate on the electron- hole 
symmetric case taking the parameters to be \x = and + U/2 = 0, where is the onsite energy 
of the interacting sites — = e^. Then the Dyson equation can be written in a A^c X Nq matrix 
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form: {d(z)} 1 = {Q Q {z)) * - E(z) with 

{^(z)}' 1 =z±-n° c - V mix (z) . 



(3.1) 



Here i is the Nq x unit matrix, and the matrices can be written in the partitioned forms; 
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S{z) 



VrF r (z) 
E u (z) E 12 {z) ... E 1N {z) 

s 21 {z) ■■. ••. ; 



(3.4) 



■ Eni{z) Unn(z). 

where l is the M x M unit matrix, and h y is the M x M hopping matrix in the y-direction; 
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(3.5) 



We note that the Hartree-Fock term of the self-energy is already included in the unperturbed 
Green's function Q (z) defined by eq. (3.1). Therefore, £{z) is the self-energy correction beyond 
the mean-field theory, which is described by the many-body perturbation theory with respect to 



We note that the mixing with the 



^c* = u Y,i=i [n^n^ - (n i < [ + n il )/2] where n ia = c\ a c ia . 
noninteracting leads is included in the unperturbed Hamiltonian, and thus S(z) depends on 
and i>r through Q°(z). In eq. (3.4) the partitioned element Uui(z) is a M x M matrix the (m, m') 
element of which corresponds to the self-energy correction between the sites located at r = (I, m) 
and r' = (l',m'), where I and m correspond to the x and y coordinates, respectively. The local 
Green's functions at the interfaces of the isolated leads, Fl(z) and -Fr(z) in eq. (3.3), depend 
on the excitation spectrum of the leads, i.e., Tt^ and Hr. We concentrate on the case i*Y = .Fr 
(= F) assuming the same excitation spectrum for the left and right leads. Then the matrices of 
the level width are given by i\(w) = —v1lmF(uj + i0 + ) and -Tr(w) = —v^lmF(uj + i0 + ). For 
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the function F(z), we consider two types models: I) semi-infinite tight-binding leads, and II) leads 
of a constant density of states. A schematic picture of the type I lead is illustrated in Fig. 3: the 
hopping matrix element is given by the nearest neighbor one, t x and t y , as that in the central 
region. Therefore the Green's function for the type I lead -F 1 satisfies a M x M matrix equation 
F\z) 



h° y -tlF\z) 



. In the case of the type II leads, we assume that the local density 
of states at the interfaces p is a constant and the bandwidth is infinity. Then the corresponding 
retarded Green's function becomes pure imaginary independent of the frequency u. Specifically, we 
consider a simple diagonal matrix of the form F (u + i0 + ) = — iirpi. Thus for the type II leads, 
the effects of the mixing are parametrized by the constant T a = irpv^ for a = L and R. 

The subband structure of the system is determined by the eigenstates of eq. fl3.5|) : h® x n = e n Xn 
for n = 1, 2, . . . , M. Due to the translational symmetry in the y-direction, the self-energy Uui(z) 
and the M x M matrix Green's function Gu/(z), which are the (1,1') partitioned element of G(z), 
can be diagonalized using the eigenstates; 

M 

(*) = £*» 4? ( 3 - 6 ) 



n=l 
M 



G w (z) = J2XnG^(z) X i. 



(3.7) 



n=l 



Therefore, the conductance eq. ( |2.9| ) can be decomposed into the sum of the contributions of the 
subbands: 

(3.8) 
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Here T^(0) is defined by r a (z) = Yl^=i Xn ^cT (0) Xn f° r a = L,R. Similarly, the Friedel sum 
rule eq. (2.11) can be rewritten as 
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(3.10) 



Here 5^ = 1 / (2i) log det corresponds to the phase shift of the n-th subband, and the 

charge displacement in each subband is given by AN$ = 2 (5 (n) /vr. As mentioned, at T = the 
conductance is determined by the value of the Green's function at u = 0, and ImX l+ (0) = 6 due 
to the Fermi-liquid property. Thus the effective Hamiltonian 7i c = 7i c + ReX' + (0), which was 
introduced in the previous section, can be used to calculate the conductance and the displacement 
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of the total charge.113' We now consider the M x M matrix T-C\ 
of Ti c . The matrix H^. w can also be diagonalized as 



CdVi 



which is the partitioned element 




M 

/ , X.n l W Xn ) 
n=l 



(3.11) 
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£n$W ~t x [8l,l'+l + +ReS}P 
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(3.12) 



where 1 < 1,1' < JV. Thus each of the modes can be mapped onto a tight-biding model in one- 
dimension with the renormalized hopping matrix element As it will be seen later, the behavior 
of eigenvalues of the N x N effective Hamiltonian defined by Ti^ ^ = {— is related to the 
transport properties. 

We calculate the conductance using the perturbation expansion with respect to U. In the electron- 
hole symmetric case, the order U 2 self-energy can be described by the diagram shown in Fig. 4; 



where 1 < < Nq- The explicit form of the unperturbed Green's function G^,(ie) can be 
obtained by taking the inverse of eq. (3.1). Since TL m \ x is included in the unperturbed part, G--, 
and Si^-, depend on the mixing matrix elements v L and f R . Note that the retarded function at uj = 
can be obtained from the Matsubara function, i.e., Et.,(0) = S^y (ie) | e _> + • The imaginary part 
of eq. (3.13) vanishes, i.e., ImS^.,(0) = owing to the Fermi-liquid property. Furthermore, due to 
the electron-hole symmetry, Re5Tt-,(0) = if j and f belong to the same sublattice. We calculate 
all the nonzero elements of Re S^,(0) carrying out the integrations in eq. (3.13) numerically. Then 
the full Green's function G^ 1 (0) is evaluated substituting the order U 2 self-energy into the Dyson 
equation. In what follows we will assume the inversion symmetry v L = t> R (= v), and take the 
parameter to be v/t x = 0.9 in the case of the type I leads. We have actually done some calculations 
also for v/t x = 0.7, but the results are similar to those for v/t x = 0.9 qualitatively. Correspondingly, 
in the case of the type II leads we set Tl = Tr ( = F ) and take the value to be T/t x = 0.75. We 
will mainly discuss the results obtained for the leads of the width M = 4. 

3.2 Noninteracting systems 

We now discuss some properties in the noninteracting case in order to show the situations we are 
considering clear. Fig. 5 shows the conductance g® d (Ep) for noninteracting electrons in the one- 
dimensional systems connected to the leads of the type I (solid line) and II (dashed line). In this 
figure, Ep is the Fermi-energy, and the size of the central region is taken to be N = 4. The coupling 




(3.13) 



8 



with the leads are taken to be v/t x = 0.9 for the type I leads, and T/t x = 0.75 for the type II 
leads. There are N (= 4) peaks of the resonance, and the peaks become sharp if the mixing matrix 
element v/t x or T/t x decreases. In two dimensional lattice with the periodic boundary condition 
in the y-direction, the conductance for noninteracting electrons is given by g® A = Sn=l 5 f id( e «) ^ n 
the electron-hole symmetric case. For instance, in the case of M = 4, the eigenvalues e n are given 
by ei = —64 = 2t y , and £2 = £3 = 0. Therefore, when t y decreases, the conductance shows peaks 
corresponding to the resonance occurring in the one-dimensional system g®Je n ). 

As an example of the the subband structure in the case of the type I leads, we show the dispersion 
relation of an ideal system E n k x = —2t x cos k x + e n in Fig. 6. Here the size along the y-direction 
is taken to be M = 4, and the hopping matrix elements are taken to be t y /t x = 1.0 (solid line) 
and t y /t x = 0.5 (dashed line). In the case of t y /t x = 1.0 the highest and lowest subbands, i.e., 
modes 1 and 4, do not contribute to the conductance because the Fermi energy is situated at the 
band edge of the subbands. These two subband become conducting for 0.0 < t y /t x < 1.0. We note 
that the modes 2 and 3 are degenerate and several curves are overlaps on the line corresponding to 
— 2t x cos k x in Fig. 6. 

3.3 Hubbard model connected to the type I leads 

In this subsection, we discuss the results for the Hubbard model connected to the type I leads. 
The conductance g for the isotropic hopping t y /t x = 1.0 is shown as a function of U in Fig. 7, 
where the size of the interacting region in the y-direction is M = 4 and that in the x-direction is 
N = 4,6, 8, and 10. The conductance decreases with increasing U. The N dependence is similar 
to that in the one-dimensional case, i.e., g decreases monotonically with increasing N showing a 
tendency toward the Mott-Hubbard insulator B As mentioned in the above, the modes 1 and 4 are 
not conducting in the case of t y /t x = 1.0, and the conductance is determined by the contributions 
of the modes 2 and 3. This situation is changed in the anisotropic cases t y /t x < 1.0. 

In Fig. 8, the conductance is plotted for several values of t y /t x , where the size of the interacting 
region is taken to be N = 4 and M = 4. There is a broad peak at finite U for 0.65 < t y /t x < 0.75. 
This shows that there are parameter regions where the total conductance increases with U even in 
the half-filled case. To clarify this behavior, in Fig. 9 the contributions of each of the conducting 
modes are plotted separately taking t y /t x to be 0.7. The resonant tunneling occurs in the modes 
1 and 4 at U/(2irt x ) — 3.3, while the contributions of the modes 2 and 3 decreases monotonically 
with increasing U. Note that the pair of the subbands whose wavenumber in the y-direction are 
k y and k y give the same contributions to the conductance. Furthermore in the electric-hole 
symmetric case, the contributions of the subbands whose eigenvalues are e n and — e n are the same. 
The occurrence of the resonant tunneling is linked with the behavior of the eigenvalues of effective 
Hamiltonian which is defined in terms of the renormalized hopping matrix element — t!j?) given by 
eq. (3.12). In Fig. 10, the eigenvalue of the N x N matrix 7i c is plotted as a function of U for the 
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mode 1 (solid line), where t y /t x = 0.7, N = 4 and M = 4. The second lowest eigenvalue crosses 
the Fermi energy at U/{2irt x ) ~ 4.6. This corresponds to the peak seen at U/{2nt x ) — 3.3 in Fig. 
9, although the values of U do not coincide. This difference is mainly due to the contribution of 
the real part of the mixing term eq. (3.3): ReVmi x (0) causes the energy shift, and the position of 
the resonance corresponding to the second eigenvalue moves toward the low-energy side. While the 

transmission probability is defined with respect to the effective Hamiltonian of the whole system 

~(n) 

, the eigenvalues of the partial Hamiltonian Tiq are useful to investigate the 
behaviors of the conductance. Note that t u , depends on the mixing matrix elements v L and u R 
because the unperturbed Green's functions used to calculate the self-energy eq. (3.13) are defined 
with respect to the whole system. In Fig. 10, the eigenvalues for the modes 2 and 3 are also plotted 
(dashed lines). In the present case, the eigenvalues of the modes 2 and 3 are the same, and the 
eigenvalues of the modes 1 and 4 are symmetric with respect to the Fermi energy u = 0. The energy 
gap between the second and third eigenvalues of the mode 2 (or 3) becomes large with increasing 
U, and it seems to show the tendency toward the Mott-Hubbard insulator. 

We have also examined the conductance in the case of M = 6, where the subbands are classified 
into two groups; e± = —€q = 2t y and €2 = €3 = —€4 = —65 = t y . In Fig. 11, the total conductance 
(solid line) and the contributions of the two groups of subbands (dashed line) are shown for t y /t x = 
0.65, where the size of the interacting region is taken to be N = 6 and M = 6. The total 
conductance show two bump like behaviors at U/(2Trt x ) ~ 1.5 and 5.1, which correspond to the 
resonance occurring in each group of the subbands. The resonance occurs at different values of U 
for the different groups of the subbands. This seems to show a general tendency that the resonance 
does not occur simultaneously in the different groups of the subbands. Thus, when the resonance 
occurs at one group, the remaining M — 2 or M — 4 subbands behave monotonically. Therefore the 
presence of the visible maximum in the U dependence of the total conductance is expected only 
for small M where the contributions of the resonating subbands are comparable to the those of the 
remaining subbands. 

The number of the eigenvalues which cross the Fermi energy is expected to increase with the size 
in the x-direction N. To clarify this feature, we calculate the conductance for N = 10 keeping the 
width M = 4 unchanged. The results are shown in Fig. 12 for several values of t y /t x . In each of 
the cases, at least one peak is visible at finite U. Specifically in the case of t y /t x = 0.5, a shoulder 
like structure is seen at U/2nt x ~ 0.5 in addition to the peak at U/2irt x ~ 3.8. To see this behavior 

precisely, the contributions of each of the subbands are plotted in Fig. 13. Two resonant peaks are 

~(n) 

seen at U/ {2irt x ) ~ 0.5 and 3.8 in the modes 1 and 4. The U dependence of the eigenvalues of 7i c 
for n = 1 (mode 1) is shown in Fig. 14. The 4th and 5th lowest eigenvalues cross the Fermi energy 
at U/(2irt x ) ~ 0.7 and 4.1, respectively. These correspond to the two resonant peaks seen in Fig. 
13 as in the case of N = 4. The energy shift caused by the mixing with the leads seems to decrease 



H qp in eq. (2.13 
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with increasing N. 



3.4 Hubbard model connected to the type II leads 

We next discuss the conductance of the Hubbard model connected to the type II leads. In this 
lead, the local density of states at the interface of the lead is a constant and the bandwidth is 
infinity. Then V m i x (u>) becomes pure imaginary and independent of u, i.e., Re V m ix(w) = and 
dV^ixi^*)/ 3°° = 0- Therefore the energy shift caused by the mixing with the leads is expected to 
be small compared to that in the case of the leads I. In Fig. 15 (a), the conductance in the case of 
ty/t x = 0.6 and 0.78 is shown as a function of U, where the size of the interacting region is chosen 
to be N = 4 and M = 4. There is a resonance peak at finite U in each line; at U/(2ntx) — 4.3 for 
ty/tx = 0.6 and U/{2nt x ) — 0.6 for t y /t x = 0.78. The contributions of each of the subbands are 
plotted separately in Fig. 15 (b). The peak seen in the total conductance is due to the resonance 

tunneling in the modes 1 and 4. Qualitatively, these results are similar to those of the type I leads. 

— ( n ) 

In Fig. 16, the eigenvalues of the effective Hamiltonian 7i c for n = 1 (mode 1) are plotted for 
(a) t y /t x = 0.6 and (b) t y /t x = 0.78. The second lowest eigenvalue for t y /t x = 0.6 becomes zero 
at U/(2irt x ) ~ 4.1. This is close to the position of the resonance peak seen in Fig. 15 (b). The 
lowest eigenvalue for t y /t x = 0.78 does not cross the Fermi energy for any values of U, although it 
is nearly zero for small U. Nevertheless, this eigenstate corresponds to the resonance peak seen at 
U/{2ixt x ) ~ 0.6 in Fig. 15 (b). This is because the lowest eigenvalue is shifted slightly up due to 
the mixing, and it becomes to cross the Fermi energy. 

We next discuss the phase shift 5^ n \ As mentioned, the charge displacement of the n-th subband 
can be deduced from the phase shift using the Friedel sum rule AiV t ^ = 2 5^ /tt, where the factor 
2 corresponds to the spin degeneracy. Specifically, in the case of the type II leads, the compensation 
theorem holds.© Then the second and the third term in the right-hand of eq. (2.10) vanish because 



V m ; x (w) is independent of uj (see also Appendix B). Thus the number of electrons in the type II 
leads is unchanged when the leads are connected to the interacting region, and the left-hand side of 
eq. (3.9), AN tot , becomes equal to the number of electrons in the interacting region ^Ziec a( c la c icr) ■ 
In Fig. 17, the phase shift is plotted as a function of U for the modes 1 and 4, where the 

parameters are taken to be N = 4, M = 4, and t y /t x = 0.78. Here the principal value of the phase 
shift is determined in the limit of U — ► by comparing with the charge displacement obtained 
independently. Reflecting the resonance seen in Fig. 15 (b), 5^ for the mode 1 starts to increase 
rapidly at U/(2Trt x ) ~ 0.6. Correspondingly in the mode 4, which is the lowest subband in this 
case, the phase shift 5^ decreases. These behaviors show how the filling of the subband changes 
when the resonance peak crosses the Fermi energy. In Fig. 17, the sum of the two phase shifts is a 
constant, i.e., 2 <JW /tt + 2 5^ /tt = 2 N with N = 4. The phase shifts for the modes 2 and 3 are the 
same, and it is a constant independent of U, i.e., 25^ /tt = N for n = 2 and 3 with N = 4. This 
means that these two subbands are half filled. Note that the total charge displacement is given by 
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AN tot = M x N in the electron-hole symmetric case. 
3. 5 Remarks 

The results presented in this section are obtained using the order U 2 self-energy given by eq. 
(3.13). Therefore, the results plotted for rather large values of U may not be sufficient in quan- 
titative sense. Nevertheless, some typical features of the results are seen even for small U. For 
instance, the resonance at finite U occurs for small value of U / (2irt x ) ^ 1 in the case of t y /t x = 0.78 
as shown in Fig. 15 (a). Therefore we believe that the qualitative features of the results, such as 
the relation between the resonance and the eigenvalues of the effective Hamiltonian, hold also for 
large U. The perturbation approach gives us a correct, in principle, description of the conductance 
and the total charge displacement with respect to the local Fermi-liquid ground state. 

§4. Summary 

We have studied the conductance through small interacting systems connected to multi-mode 
leads based on a local Fermi-liquid approach. At T = 0, the conductance and the total charge dis- 
placement are determined by the value of the Green's function at the Fermi energy. Since the exci- 
tations at the Fermi energy do not decay at T = 0, there exists the one-particle Hamiltonian which 
reproduces these two quantities exactly. We have applied this formulation to a two-dimensional 
Hubbard model of finite size in the electron-hole symmetric case. We have calculated all the matrix 
elements of the order U 2 self-energy and determine the parameters of the effective Hamiltonian up 
to terms of order U 2 . Specifically, we have examined two types of the noninteracting leads: I) semi- 
infinite tight-binding leads, and II) leads of a constant density of states. The results are similar, 
qualitatively, in both types of the leads. The conductance shows maximums in the U dependence 
for some ranges of t y /t x , where t x and t y are the hopping matrix element in the x- and y-directions, 
respectively. This means that there exists parameter regions, where the total conductance increases 
with U, even in the half-filled case. By decomposing the total conductance into the sum of the 
contributions of the subbands, it is clarified that the peaks of the conductance are caused by the 
resonance occurring in some group of the subbands which have the similar symmetric properties. 
The phase shift of the subbands obtained from the Friedel sum rule shows a typical changes when 
the resonance occurs. The resonance generally does not occur simultaneously in different groups of 
the subbands, and the subbands of off-resonance behave monotonically. Therefore the maximum 
in the U dependence of the total conductance is expected only for the mesoscopic systems in which 
the number of the conducting modes is small enough and the contributions of the subbands of 
on-resonance are comparable to the those of the remaining subbands. 

The perturbation approach examined in this work can be applied to interacting electrons in 
disordered systems. Especially, the analysis of the eigen values of the effective Hamiltonian, which 
is used in the present study to investigate the behaviors of the resonance, may be combined with 
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the Thouless-numbeinH) and finite-size scalingEiP methods. The application along this line seems to 
be interesting in relation to the metal-insulator transition observed in two-dimensional systems® 
Furthermore, extensions to the finite temperatures© and nonequilibrium situations!^ are left for 
future studies. 
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Appendix A: Conductance 



In this appendix, we provide the derivation of the dc conductance in the multi-mode case eq. (2.9) 



by generalizing the derivation given in the single-mode case© The current operator corresponding 
to the mixing term eq. (2.6) is given by 

m l 



J L = ie EE U L, m (4 mCT C £ m a- C La C X mCT ) . ( A4 ) 

Mr, 

J R = 1 e E E %,m ( 4i m a c N m * - c N m * c n m a ) • (A-2) 

777=1 a 

Here Jl is the total current flowing into the sample from the left lead, and Jr is the current 
flowing out to the right lead from the sample. These currents and the total charge in the sample 
Pq = e J2jec,a c jcr c ja satisfy the equation of continuity dp c /dt + Jr. — Jl = 0. In the linear response 
theory, the dc conductance is given by 

K aa ,(m) = [ dr ( T T J Q (t) J a i (0) ) e i UlT , (A-4) 
Jo 

where a, a' = L,R, and K~^ a ,{oj) is the retarded function which is obtained from the Matsubara 
function K aa i(ivi) by the analytic continuation if/ — > to + i0 + . The dc conductance eq. (A-3) 
corresponds to the w-linear imaginary part of K7~ ,(uj), and it is independent of a and a' owing to 
the current conservation.!!© Moreover, K Q / a (z) = K Qa /(z) because of the time reversal symmetry, 
and thus the imaginary part of , (to) corresponds to the discontinuity of K aa > (z) along the real 
axis in the complex z-plane. Therefore, the dc conductance is equal to the coefficient of the v sgn v 
term of K aa >(w where sgn^ is the sign function. In what follows, we extract this singular term 
from K aa i(iv) taking the current operators to be a = R and a' = L. 

At T = 0, -Krx^) is written, treating the Matsubara frequencies to be continuous variables, as 
[see also Fig. 1] 
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4lV) 



-°° de 



(A-5) 



xTr 



/+oo ^ £ 
-oo 2^ 

^ R Gnc(te + iv) v L GxAf{ie) + t> R G^j(ie + iz/) w L G^(ie) 



-v R G n i(ie + iz/) t> L Gcjv(is) - v R G^ c {ie + iz/) t> L G m (ie) 
<7<t' {j}ec ' 

x | G C £(ie + iv) v L G IC (ie) - G C i{ie + iz/) v L G C cfe) > 

{ ' 3134 



(A-6) 



x \ GcAr(ie') v R Gnciye' + iv) - G C n{^') v R G^die' + iz/) 



(A-7) 



3332 



Here i> L and u R are diagonal matrices corresponding to -u L m and t> R m , respectively. The subscript 
C denotes a set consisting of Nq sites in the central region. Gcc(z) is a Nq x Ml matrix, and the 
(j,m) element is given by Gj : c m (z) with j G C. Also, Gkc(z) is a M R x Afc matrix, and the (m,j) 
element is given by Gfi m j(z) with j G C. The matrix Green's functions in eqs. (A-6) and (A-7) 
are defined in this way. r^, 2 '^ 4 (i£i, 1E2; 163, 1^4) is the vertex corrections due to the inter-electron 
interaction, and illustrated in Fig. 18. Since we are considering the interaction which is switched 
on only in the central region, the Green's function satisfies following relations at the interfaces; 



G-ji-yiz) = - F R (z) v R G^{z) for 7 = £, 1, C, Af 
G lC (z) = - G lX {z) v L F L (z) for 7 = J, C, jV, K 

Using these relations, eqs. (A-6) and (A-7) are rewritten as 



v R [ F R (ie + iz/) - F R (ie) ] v R Gxj(ie + iz/) 



x v L [F h (ie + iz/) - F h (ie) ] v h G IAf (ie) 



x J Gaiie + iz/) v L [F L (k + iz/) - F L (k) ] v L G IC {ie) 



(A-9) 
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x I GcAr(ie') v R [ F R (ie' + iv) - F R (ie') ] v R G Mc {\e' + . (A- 10) 

We can now extract the singular vsgnv term of Arx^) using eqs. (A-9) and (A-10), as in the 
case of the single Anderson impurity]!! El* It can be shown that there is no singular v sgn v term in 
K R b l(iv) due the two factors [F^(ie + iu) — F^(ie)] and [F^(ie' + iv) — F^(ie')] which have different 
frequencies. Thus, the vsgnv term comes only from K$(iv). It can be obtained from eq. (A-9) 
and yields eq. (^9|). 

Appendix B: Priedel sum rule 

We provide here the outline of the derivation of the Friedel sum rule in the single-mode case eq. 
(B-7) following Langer and Ambegaokar.© The Green's function in the left lead ij G L is written 



as 



G+M = Ft d (u) + F^) v L G+(u) v L F^(uj) , (B-l) 

F ^ ) -\u J -^ + K + - (B ' 2) 

Here Fjj(uj) is the retarded Green's function of the isolated lead, and e n ^ and 0„l(z) are the 
eigenvalue and the corresponding eigenfunction of the one-particle Hamiltonian TL-^. Using the 
orthogonality relation J2ieL 4> n lW ^n' L W = ^rm'i we obtain 

dFk 



M it 



iGL 



00 

duj 



v h Gl lVh . (B-3) 



Here Fq = F£, which is the Green's function at the interface. Using eq. (B-3) and the correspond- 
ing relation for the right lead, the displacement of the total charge defined by eq. (2.10) is written, 
at T = 0, as 



AN, 



tot 



-Im f° da; I £G+ + £[G*-i$] + 



7T 



R 



2 Z" 
— Im / da; Tr 

7T 



_ + mix 

9a; 



Im / da; ( — — Tr 

vr J-oo \ ooj 



log{^ + }^ 



Tr 



(B-4) 



Here we have used a matrix notation as that was used in the text eqs. (3.1)-(3.4) for the multi- 
mode case: the Green's function in the central region is denoted by Q(z) = {Gij(z)} with ij G C and 
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the Dyson equation is given by {G(z)}~ 1 = zx — Hq — V m i x (z) — £(z), where H.q = {— 1§ — fiSij}, 
and Vraix(z) and £(z) are the single-mode version of eqs. (3.3) and (3.4), respectively. 

It has been shown that the contribution of the second term in the third line of eq. (B-4) van- 
ishes.0® Thus, eq. (B-4) can be rewritten as 

L log[det{g^(0)}" 1 /det{g + (0)}" 1 " 



AN, 



tot 



7T1 
1 

7ri 



log det [i + 2ilm{v+ x (0)} g + (0) 



(B-5) 
(B-6) 



For obtaining eq. (B-6), we have used the Fermi-liquid property Im-E^O). The matrix V m ix has 
only two non-zero elements as eq. (3.3). Therefore, the determinant in eq. (B-6) can be written in 
terms of a 2 x 2 matrix as 



AN, 



1 



tot 



7T1 



log [ det S] , 



(B-7) 





"1 


0" 




s = 






- 2i 







1_ 





r L (o) o 
o r R (o) 



G+(0) G+ N (0) 
GnM G% n (0) 
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Fig. 1. Schematic picture of the system. 
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Fig. 2. Schematic picture of a 2D Hubbard model: (•) interacting sites. 
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Left lead 

Fig. 3. Schematic picture of the type I lead. 
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Fig. 4. The order U 2 self-energy Ejy(ie). 
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Fig. 5. Conductance vs. Ef for noninteracting electrons on the one-dimensional chain of the size N — 4 connected 
to the leads; (solid line) type I leads and (dashed line) type II leads. 
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-71 71 

Fig. 6. Subband structure of the type I lead of the width M = 4, where (solid line) t y /t x = 1.0 and (dashed line) 

I t x = 0.5. 




U/(2tc1 x ) 

Fig. 7. Conductance vs. (7 in the isotropic case of t y /t x = 1.0. The size of the interacting region is chosen to be 
M = 4 and N = 4, 6, 8, 10. Type I leads. 
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U/(2tc1 x ) 

Fig. 8. Conductance vs. U in the anisotropic case t y /t x < 1.0. The size of the interacting region is chosen to be 
M = 4 and iV = 4. Type I leads. 



22 




1 2 3 4 5 6 

U/(2tc1 x ) 



Fig. 9. Conductance of each subband, where M = 4, TV = 4, and ^/£ x = 0.7. Type I leads. 
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Fig. 10. Eigenvalues ol Tic for n = 1 (solid line) and n = 2, 3 (dashed line), where M — 4, N = 4, and t^/t^ = 0.7. 
The origin of the energy is set to be Ef = 0. Type I leads. 
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Fig. 11. Conductance vs. (7 for the system of the size M — 6 and iV = 6: (solid line) total conductance and (dashed 
line) contributions of each mode for t y /t x = 0.65. Type I leads. 
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Fig. 12. Conductance vs. £7 for the system of the size M — 4 and N = 10, where t y /^ = 0.4, 0.5, 0.6, and 0.7. Type 
I leads. 
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Fig. 13. Conductance of each subband for the system of the size M = 4 and N — 10, where t y /t x = 0.5. Type I 
lead. 
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Fig. 15. Conductance vs. U for the system of M — 4 and N = 4 in the case of the type II leads: (a) total conductance, 
and (b) contributions of each mode. The hopping matrix element is taken to be (dashed line) t y /t x = 0.6, and 
(solid line) 0.78. 
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Fig. 16. Eigenvalues of Tic f° r n = 1, where M — 4 and N = 4. The hopping matrix element is taken to be (a) 
tyjt x = 0.6 and (b) 0.78. The origin of the energy is set to be Ef = 0. Type II leads. 
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Fig. 18. The vertex function r^, 2 '^ 4 (ie + i^, ie' + i^; ie', ie). 
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